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Abstract 

Laboratory experiments point out the existence of patterns made of alternately 
laminar and turbulent oblique bands in plane Couette flow in its way to/from turbu- 
lence as the Reynolds number R is varied. Many previous theoretical and numerical 
works on the problem have considered small aspect-ratio systems subjected to pe- 
riodic boundary conditions, while experiments correspond to the opposite limit of 
large aspect-ratio. Here, by means of fully resolved direct numerical simulations 
of the Navier-Stokes equations at decreasing i?, we scrutinize the transition from 
temporal to spatiotemporal behavior in systems of intermediate sizes. We show that 
there exists a streamwise crossover size of order Lx ^70-80 h (where 2h is the gap 
between the plates driving the flow) beyond which the transition to/from turbulence 
in PCF is undoubtedly a spatiotemporal process, with typical scenario '[turbulent 
flow] [riddled regime] [oblique pattern] [laminar flow]', whereas below that 
size it is more a temporal process describable in terms of finite-dimensional dynam- 
ical systems with scenario '[chaotic flow] [laminar flow]' (via chaotic transients). 
In the crossover region, the '[oblique pattern]' stage is skipped, which leads us to 
suggest that an appropriate rendering of the patterns observed in experiments needs 
a faithful account of streamwise correlations at scales at least of the order of that 
crossover size. 



1 Introduction 

During the last few years there has been a resurgence of interest regarding the formation 
of laminar-turbulent patterns in wall-bounded shear flows having very large aspect ratios 
(when the lateral dimensions, along x and z, are more than two orders of magnitude larger 
than the relevant wall-normal dimension, along y). The phenomenon was discovered in 
circular Couette flow (CCF) by Coles and Van Atta j8l[9] who called it 'spiral turbulence'. 
The pattern can be observed when the two cylinders of CCF rotate in opposite directions 
in a specific velocity range, which is only a part of the complete bifurcation diagram 
[1]. Only one or two stripes of the laminar-turbulent pattern were observed in those 
experiments owing to their low aspect ratio (the ratio of the perimeter to the gap between 
the cylinders). It is only later that Prigent et ai [32j made detailed measurements of the 



1 



Figure 1: Top: Schematic of PCF with top and bottom plate moving at ±[/p (at y = 
±/i) separated by a distance 2h. Bottom: Experimental visualization of laminar (dark)- 
turbulent (light) pattern in PCF at i? = 358, = 770 h, = 430 h (courtesy A. Prigent). 
For comparison, most of S. Bottin's experiments were performed with = 284: = 
72.6 h [HI [5]. The size of the MFU is A^^ = 27r/a, = 27r/7, where a and 7 are the 
fundamental wavenumbers in the streamwise and spanwise directions, respectively. Here 
we take the values corresponding to the self-sustained exact solution found by Waleffe 
[IDl, a = 0.49 and 7 = 1.5 and thus A^ = 12.8 and A^ = 4.2. 

relevant features of this phenomenon in their very-large aspect ratio setup, where about 
10-15 stripes could be obtained. 

A similar event of laminar-turbulent pattern formation was also observed by the same 
team [32] in large- aspect-ratio plane Couette flow (PCF), a schematic of which is shown 
in the top panel of flguredj the bottom one displaying a snapshot of the pattern obtained 
experimentally. When the control parameters of CCF and PCF are made equivalent 
by using appropriate scales, both CCF and PCF display patterns in almost the same 
parameter range [22j, with the obvious difference of streamwise periodicity present in 
CCF and not in PCF. In the following we leave aside the case of CCF whose bifurcation 
diagram is slightly more complicated owing to the interplay of centrifugal instability 
mechanisms and concentrate our attention on PCF in the range of parameters relevant 
to the transition to turbulence. 

Traditionally, the control parameter is the Reynolds number, and for PCF it is deflned 
as R := Uph/jy^ here based on the half channel height /i, the plate velocity [/p, and the 
kinematic viscosity in the following h and the advection time Ta = h/Up are used as 
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length and time units, respectively. (The viscous relaxation time is Ty = h'^/iy^ so that 
R = Ty/r^.) The transition process can be studied by increasing or decreasing R. The 
essential point is that the laminar flow is stable against inflnitesimal perturbations for all 
i?, so that a direct transition to turbulence is observed when flnite amplitude perturbations 
are introduced, according to a globally subcritical scenario (for an introductory review 
see jSU Chap.7]). Studies at increasing R are strongly sensitive to the amplitude and 
shape of the triggering perturbation. Quantitatively diflFerent results about the transition 
can accordingly be obtained under diflFerent protocols, e.g. [38[ [5], while the picture of 
course remains qualitatively unchanged. In contrast, the experiments of Prigent et al. (see 
[3T] for details) , were systematically performed by varying R in small steps while waiting 
for statistical equilibrium at each i?, which helps us to clearly identify several stages. 
First, beyond Rt ^ 410, turbulence is essentially uniform or 'featureless', borrowing the 
term introduced by Andereck et al for CCF [Ij. Next, oblique laminar-turbulent bands 
appear upon decreasing R slowly below Rt. The amplitude of the laminar-turbulent 
modulation grows continuously as R is further decreased, which has been interpreted as 
a supercritical bifurcation in the presence of a noise reminiscent of featureless turbulence 
[32j| . The bands next become fragmented and turn into irregular oblique turbulent patches 
which seem sustained for R> R^ ^ 325 but decay in a flnite time for R < R^. The value 
of i?g mentioned above has also been obtained in experiments where turbulent spots were 
triggered P, ^ . Below i?g the lifetimes of turbulent spots are distributed according to 
decreasing exponentials whose characteristic time is seen to diverge as R approaches Rg 
from below [6j. Finally, for R < R^ ^ 280 large perturbations relax without measurable 
waiting time and in a mostly monotonic way. 

Direct numerical simulation (DNS) of these large aspect ratio systems was delayed 
because of the huge computational requirements. Barkley and Tuckerman [2j were the flrst 
to obtain the band patterns in fully resolved simulations of the Navier-Stokes equations, 
however in carefully chosen narrow and tilted computational domains. Their choice of 
domain however precluded the occurrence of patterns with defects or orientation changes 
inside the flow. This restriction was overcome in the DNS by Duguet et al. [13j who 
recovered the experimental flndings of Prigent et al. in a fully resolved very large aspect 
ratio system. Similarly, the spiral regime of CCF was numerically obtained by Meseguer et 
al. [27] and Dong [llj. An oblique band pattern was also found numerically by Tsukahara 
et al. [39] in the case of plane Poiseuille flow (PPF). The above three systems, CCF, PCF 
and PPF are prototypical for the study of laminar-turbulent patterns; all of them are 
conflned by two walls in direction whereas the other two dimensions extend to inflnity, 
in the limit. 

Two other flows of theoretical and technological importance showing behavior similar 
to the above-mentioned systems are pipe flow and the flat-plate boundary-layer flow. Pipe 
flow is in some sense simpler, owing to only one spatially-extended dimension, the length of 
the pipe. Inside a speciflc Reynolds number range, localized patches of turbulence ('puffs') 
are observed separated by laminar regions. At present it is not clear how they are related 
to the oblique bands in other systems [28j. On the other hand, the boundary layer flow 
case gets more complicated due to the fact that the relevant wall normal dimension, the 
boundary-layer thickness, increases with the flow direction and, with it, the relevant value 
of R. Turbulent spots that appear amid laminar flow during the transition are however 
strikingly similar to those growing into oblique bands observed in PCF. 

In contrast with instabilities in closed flow systems such as Rayleigh^Benard convec- 
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tion for which the aspect ratio is a genuine control parameter that can be varied from 
smaU to large by changing the position of walls in all directions [241 Chap. 3], open flows 
develop most often in domains physically less constrained by lateral walls, while stream- 
wise boundary conditions make it reasonable to accept the assumption of translational 
invariance in that direction, at least locally. In this respect, all of the investigations men- 
tioned above, both experimental and numerical, are related to large aspect ratio systems 
and represent the typical situation, even if the experimental setup is difficult to construct 
|38[ |5[ |32] or the DNS computationally demanding [131 EH]. Despite its relevance, this 
is not the situation that has been considered in many recent works on the transitional 
problem. Indeed, most of the theoretical and numerical work has been performed in con- 
figurations confined by periodic boundary conditions at small distances. Seminal work 
in the 1990's bore on high resolution simulations devoted to the identification of the so- 
called minimal flow unit (MFU) below which turbulence can not be sustained, as done by 
Jimenez and Moin [18], later used for the elucidation of the self-sustaining process (SSP), 
the mechanism of sustenance of turbulence by which streamwise vortices induce streaks 
that break down to regenerate the vortices, following Hamilton et al. [17j, next to the 
discovery of exact coherent states of Navier-Stokes equations, following Nagata [29| and 
Waleffe [40j , and to the study of the boundary of the attraction basin of the laminar fiow 

In fact setting boundary conditions at small distances reduce the infinite-dimensional 
dynamical problem posed by the Navier-Stokes equations to a finite-dimensional problem 
as long as R remains moderate, which is the case in the transitional regime (but would 
be insufficient in the high-i? limit where fine vortical structures develop already at the 
level of the MFU). Accordingly all the works in small domains to some extent have come 
under the purview of temporal dynamics and chaos theory. The investigation of such 
small systems has indeed been extremely valuable [T^ [T6] . with important achievements 
such as the recent findings of 'edge states' and 'localized solutions' in PCF [121 EI], as 
an encouragement to make the connection to special solutions and associated bifurcation 
structures obtained in pattern- forming model equations [71 [36]. 

However, any accurate representation of the dynamics at the level of the MFU, though 
remaining instructive, is not informative of the experimental situation since the smallest 
setups that have been used should rather be analyzed as two-dimensional arrays made of 
tens or hundreds of MFUs (see caption of Fig. [1]), which allows for global spatiotemporal 
dynamics, while placing periodic boundary conditions at the scale of the MFU grossly 
overestimates the coherence of the fiow. The size of the system MFU) and the co- 
existence of two possible local states, either laminar or chaotic, each corresponding to 
a possible temporal regime at the MFU scale, make it possible for whole regions, either 
laminar or turbulent, to coexists in physical space. The possibility of such modulations 
roots the recourse to concepts from the theory of spatiotemporal chaos [lOj. Within this 
general framework, the pattern featureless transition in PCF, TCF and other similar 
flows, appears to be a symmetry-restoring bifurcation observed upon increasing R for 
which an order parameter can be defined ^21 H], with the understanding that the base 
state is the translationally invariant, strongly noisy, featureless regime beyond R^. 

A previous attempt to reach the spatiotemporal level directly via modeling by one of 
us [20j, though promising [23j, however failed to reproduce the bands, due to insufficient 
wall-normal resolution. Another avenue to spatiotemporal dynamics is through DNS of 
the Navier-Stokes equations but, as already mentioned, long duration simulations of wide 
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Figure 2: (Color online) Schematic of the expected bifurcation diagram for PCF including 
the locations of domains with the size of the order of MFU and that of pattern-forming 
systems. 



enough domains is still too demanding. Reducing the wall-normal resolution has been 
shown a viable option at a qualitative level \26l but it was not clear that the quantitative 
shift observed on the transitional range [i?g,i?t] was without hidden consequence on the 
pattern formation problem. 

The aim of the present work is at studying, by means of fully resolved DNS of Navier- 
Stokes equations, what happens when the system size L (to be defined more precisely 
later) increases and the temporal dynamics gives way to a spatiotemporal one. The 
next section anticipates the outcome of the study before giving details on the numerics. 
Section [3] contains the main results extracted from the numerical simulations as well as 
the bifurcation diagram. The final section summarizes the study and draws conclusions. 
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2 General framework for the study of transitional 
PCF 

2.1 Expected bifurcation diagram for the Turbulent-Laminar 
transition 

Figure [2] gives a schematic three-dimensional view of the bifurcation diagram for tran- 
sitional PCF in the two well-studied cases of MFU-like systems and large aspect-ratio 
systems, and propose an interpolation between these two cases as an educated guess. We 
use i?, and L as a coordinate system. Besides i?, the Reynolds number, is a measure 
of the distance to laminar flow, e.g. the time-averaged, domain- averaged rms of the ve- 
locity departure from the laminar flow proflle, and L is the characteristic lateral domain 
size. The diagram shows two thick black lines in two (R-D) planes, one at L = (9(MFU), 
the other at L = (9(100 MFU) where '100' is just meant to be a large number typical of 
laboratory experiments. They sketch the variation of D as a function of i?, as the latter 
parameter is gradually reduced from a high value where the system is turbulent. 

The flrst curve, at L = (9(MFU) is typical of a subcritical bifurcation with an abrupt 
jump from turbulence (D flnite and large) to laminar flow (D = 0). Here it is just 
a tentative sketch of the ideal situation where temporal chaos would break down in a 
single stage via attractor crisis [t30j. This guess is consistent with the observations of 
exponentially distributed lifetimes of chaotic transients and the divergence of the average 
lifetime as the putative crisis point is approached [H]. In fact, the very existence of a 
crisis at flnite R (here situated around R = 410) has not been proved so far but this does 
not change our picture drastically. 

In contrast, the second line at L = (9(100 MFU) corresponding to large aspect ratio 
systems is supported by laboratory and computer experiments [32^ 13j. As R is lowered 
below i?t ^ 410, the distance D begins to decrease owing to the coexistence of laminar 
and turbulent regions in variable amounts in the whole domain. Whereas the turbulence 
intensity seems to decrease only slightly with i?, the main part of the decrease of D, 
a quantity averaged over the whole surface of the system, has to be attributed to the 
increase of the laminar fraction, the complementary of the turbulent fraction. The spatial 
organization of this laminar-turbulent coexistence is expressed by the term 'pattern' 
(henceforth shortened as T') used in the flgure. At i?g ^ 325, the distance D drops to 
zero because the regime observed in the long-time limit is laminar, but this does not mean 
that one cannot observe transiently turbulent patches with smaller turbulent fraction and 
exponentially distributed lifetimes [6j. 

Thick red dashed lines in Fig. [2] suggest changes in the diagram as L varies, showing 
up two regions more: 'LAMINAR' ('L') and 'turbulent' ('T') or 'chaotic' ('C'). When 
L is large, 'T' actually means 'featureless turbulence' but, when L is small, 'temporal 
chaos' would better reflect the spatial coherence in the flow, hence the 'C. This coherence 
forbids the emergence of sub-domains that could be identifled as laminar or turbulent in 
the system, which is no longer the case when L gets larger. For reasons of topological 
continuity, there should be a crossover size (shown by a hatched band) below which PCF 
displays two states, 'L' and 'C, and above which three states, 'L', 'P', and 'T', the 
hatched band extending into the 'C/T' region to mark the change from temporal chaos to 
spatiotemporal chaos. It will shown below that, such a bifurcation diagram indeed exists. 
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Table 1: Computational domains and the corresponding grid points 
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with however some peculiarities in the T' region. 
2.2 Numerical simulation details 

Direct Numerical Simulation (DNS) of the Navier-Stokes equations is carried out us- 
ing Gibson's well tested, freely available DNS software ChannelFlow^ [15j. It is a 
pseudo-spectral code using Chebyshev polynomials in the wall-normal direction with no- 
slip conditions at = ±1 and in-plane Fourier modes adapted to periodic boundary 
conditions at distances L^^^ in the (x, z) directions. The number of Chebyshev modes 
is A^^ = 33, which is well suited to resolve all the relevant modes of turbulent flow in 
the range of R studied according to Duguet et al [13] (see also note 0) . The time step 
is such as to keep the CFL number between 0.4 and 0.6. Some details of the various 
domains used in simulations are given in table [U where A^^^ and A^^ are the number of 
collocation points in x and z directions, respectively, and and L^, the corresponding 
domain lengths. Since the 3/2 rule is applied in all cases to remove aliasing, this cor- 
responds to solutions in the Fourier space using N'^ ^ = \^x,z modes, or equivalently to 
space steps L^^z/N'^^ = ^Lx^z/]^x,z- Domains are chosen with increasing size, measured 
here by the diagonal length L := + , which will also be used in referring to various 
simulations. Moreover, in all cases 9 := tan~-^ (L^/L^^) is kept between 20-30°, consistent 
with previous investigations [32", ^3]. For certain domain sizes, extra computations with an 
increased number of grid points are also carried out to ensure the independence on grid 
resolution, e.g. for L = 16.9, A^^^ was raised to 96, and for L = 65.4, (A^^^, A^^) were raised 
to (282,128) without any signiflcant quantitative changes in the results. 

The laminar base flow simply reads Viam = where x is the unit vector in the 
streamwise direction. The velocity perturbation v is obtained as v = v — t/x, where 
V := + vy + w z is the instantaneous velocity fleld solution to the Navier-Stokes 
equations that are simulated. For further reference. Figure [3] displays typical snapshots of 
the solutions obtained for R = 410 in the different domains that we consider. We illustrate 

^This corresponds to an average wall- normal spacing Ay~^ = 1.81 based on Rr = 32; the superscript 
denotes quantities scaled by the viscous length unit v/ur, where Ur is the friction velocity ^Jryjj p ^ 
Tyj being the shear stress at the wall and p the fluid density; further, is defined as Urh/v and R^ = 32 
roughly corresponds io R = 420. 
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Figure 3: (Color online) Streamwise velocity fluctuations {u) in the plane y = for 
different domains at i? = 410: L = 143.1 (top-left), 109.6 (top-right), 98.5, 87.3, 76.2, 
65.4, 35.3, 25.6, 16.9 (bottom). 
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the fluctuation fleld using u evaluated in the plane ?/ = 0, which is a direct trace of the 
perturbation since viam = there. The different domains are displayed, starting with 
L = 143.1 and decreasing the size from top-left to right-bottom. The snapshots presented 
are approximately proportional to the domain sizes. Colors red and blue correspond to 
u ^ ±0.5. These images show that the very large structures discovered by Komminaho 
et al [19], in a domain of size 287r xSn {L 91.5) at i? = 750 deep inside the featureless 
regime are already present at i? = 410, in the immediate vicinity of Rt. In all cases the 
flow flelds are dominated by streamwise elongated (or streaky) structures which are more 
or less alternating in the spanwise direction. In the smallest domains, the development 
of these streaky structures is severely constrained. In larger domains, i.e. for L = 76.2 
and larger, pockets of laminar flow can be observed already at i? = 410. More details 
about these structures and their signiflcance as R is decreased will be discussed in the 
subsequent sections. 

The results to be described now, in connection to the temporal/spatiotemporal issue 
in domains of sizes varying from a few units of MFU to ones where patterns or bands 
appear, are all obtained using an 'adiabatic protocol' to be described below according to 
which, starting from a turbulent state at high Reynolds number, R is reduced by steps of 
AR and the simulation is run for AT, repeatedly down to the laminar regime. 



3 Results 

3.1 Fluctuations in varying domain size 

Here we use two global measures of fluctuation intensity, the overall rate of energy dissi- 
pation per unit volume: 

V := j [ ^ [ [ ^ {\\/u\^ + \\/v\^ + \\/w\^) dxdydz, (1) 

LxLyLz Jq J_l Jq 



where = 2, and the rate of mechanical energy input: 




LxLyLz 



du 

y=i 



dx dz . (2) 



They are normalized such that, for the laminar solution, both V and I are equal to 1. For 
any flow fleld, on average V is equal to X, and their instantaneous diflFerence is a measure 
of the fluctuation of the system's energy: 

£{t):=-^— r I [ ^ lMx,y,z;t)\fdxdydz, (3) 

J^xJ^yJ^z Jo J -I Jo 

SO that, directly from the Navier-Stokes equations, one derives d£/dt = R~^{I — V). 

Figure m shows the trajectories followed by the system as projected on the V-I plane 
for four increasing domain sizes L: 16.9, 35.3, 76.2, and 143.1. An adiabatic protocol is 
followed, similar to the annealing experiment performed by Schmiegel & Eckhardt [3^ 
in a MFU-sized system or Barkley & Tuckerman [2j in their oblique domain. We start 
at i? = 450 for the smallest domain and at i? = 420 for the others. By adiabatic, 
we mean that R is reduced by steps AR every AT, the flnal state at a given R being 
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Figure 4: (Color online) Input X vs. dissipation P as a function of time at different R; 
L = 16.9 (a), 35.3 (b), 76.2 (c), and 143.1 (d). See text for details. 
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Figure 5: (Color online) Distance to laminar flow (t/rms) at diflFerent R for L = 16.9. 
Adiabatic protocol with time series of duration AT = 5 • 10^. 

used as initial condition for the simulation at i? — AR. Usually we take AT = 10^ and 
AR = 10 for R > 320 and AR = 5 below 320. Not all the traces are shown in order 
not to over-print the flgure. Traces tending to the (1,1) point in the V-I plane express 
'turbulent laminar' breakdown, which is the case for R = 400 when L = 16.9, R = 370 
for L = 35.3, etc. This is also the case for R = 300 when L = 143.1 but the simulation has 
been interrupted before complete relaxation toward the laminar regime. For the smallest 
domain, it is observed that the trajectory quickly falls in the neighborhood of the diagonal 
(V = I) and then continues to revolve in a region along this line. Decreasing R results 
in a slight decrease in the wandering but the point representing the state of the system 
stays along the diagonal and close to it. As the domain size is increased, the amplitudes 
of the excursions along the diagonal, and away from it, decrease and the dynamics seems 
slower; the transition from one R to another also seems smoother. These characteristics 
can be understood as the result of averaging over larger domains, in connection with the 
extensive character of the featureless turbulent regime examined in §3.21 below. 

A complementary piece of information is obtained from the distance to laminar flow ^ 
whose time average was proposed in Fig.[2]to characterize the flow regime on a global scale. 
As a relevant measure of the distance, we take the volume- averaged root-mean-square 
value of the velocity perturbation v, here denoted t/rms- Time series of this quantity 
during simulations performed according to the adiabatic protocol described earlier are 
shown in Figs. [5] and [7] in two representative cases, L = 16.9 and L = 143.1, respectively: 

• As seen in Fig. [5]for L = 16.9, while the average of t/rms for diflFerent R does not 
vary much, at some reduced R the flow suddenly relaminarizes. This process was found 
to be probabilistic, e.g. by Schneider et al ^35j, and explained as arising from chaotic 
transients associated to the dynamics around a tangle in phase space. Snapshots taken 
during the decay of the chaotic state at i? = 380 are displayed in Fig. El bottom. Whereas 
images taken before decay look similar to the one at t = 2000, breakdown of the chaotic 
state is seen as a fast damping of small scale structures {t = 2050), leaving just a pair 
of streaks {t = 2100, 2150) that progressively fade away {t = 2250, 2350). Here decay is 
observed for R = 380 whereas a similar decay happened at i? = 400 in the experiment 
reported in Fig. H^, with diflFerent initial conditions and a different protocol (AT = 10^ 
instead of 5 • 10^ here) but this merely expresses the probabilistic character of chaotic 
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Figure 6: (Color online) Top: Distance to laminar flow (t/rms) during the end of the 
transient at i? = 380 for L = 16.9. Bottom: Snapshots during decay; time, corresponding 
to points on the graph, is indicated below the images; exceptionally, the streamwise 
direction is along the vertical. 
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Figure 7: (Color online) Distance to laminar flow (t/rms) at different R for L = 143.1 
{9 = 26.6°). Inserted flgures show snapshots of u at various R (indicated on top-left) 
after simulation over AT = 10^. 

transients whose lifetimes are strongly sensitive to initial conditions [14j. It should be 
noted that the exact size of the domain certainly matters for such small systems since, 
using their annealing protocol at a velocity corresponding to the rate AR/ AT = 2 x 10~^ 
of Fig. El Schmiegel & Eckhardt found that transient chaotic dynamics with lifetime 
of order 2500 time units could be maintained down to R 310 in a system of size 27r x tt, 
whereas in a system of size 57r x 27r we obtain short-lived chaotic transients already at 
R = 380 (though we have not performed a detailed statistical study). On the other hand, 
a deflnite mark of spatiotemporal dynamics will bear a much weaker sensitivity to the 
precise in-plane dimensions of the system. 

• When L = 143.1, upon further lowering of R below 400, the average value of t/rms 
is seen to decrease regularly with R until relaminarization occurs, here for R = 300. Also 
shown in Fig. [71 are the snapshots of u at various R (indicated on top-left) at the end of 
a simulation step (AT = 1000). At R = 380, laminar troughs^ i.e. small patches where 
turbulence seems depleted can be observed (and not yet patterns). Upon further reduction 
of R the troughs join to form a well deflned pattern of oblique bands, alternatively laminar 
and turbulent (snapshot at i? = 360). The subsequent decrease in [/rms is due to the 
increase in the width of the laminar band. The flnal decay to laminar flow is neither 
sudden nor uniform in space as in small domains, but, through a local breaking of the band 
leaving separated turbulent patches (see snapshot at i? = 300) that recede and disappear. 
Only at the very last stage, when the turbulent patches have reached a small size, they do 
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collapse under the effects of viscous dissipation like the flow in small domains. All along 
the decay, the internal structure of the turbulent spots is not fundamentally different from 
that of spots growing from a flnite-amplitude localized disturbance at higher Reynolds 
numbers, as obtained experimentally [381 El [6] or numerically [2Tl[13]. It should be noted 
that, here, decay is observed at i? = 300 while experiments at large aspect ratio would 
rather suggest decay for R < ^ 325. A flrst reason for this delayed breakdown may be 
that periodic boundary conditions tend to stabilize the pattern. A second, certainly more 
important, reason is that according to our adiabatic protocol, R is changed by AR every 
AT and that AR may be too large and AT too small. Basically, the systems 'knows' 
at which value of R the simulation is running only after several viscous relaxation times 
Ty = Rts,^ where Ta is the advection time /i/t/p, which is also our time unit. Accordingly, 
AT = 10^ ^ 3rv, which is barely sufficient to reach the steady state. This is all the 
more true that, owing to the subcritical character of the laminar-turbulent coexistence, 
we should allow for long waiting times implied by the nucleation of the stable (laminar) 
state within the metastable (turbulent) state at the origin of the breaking of the band 
that causes the decay [25] . The need for long lasting simulations in wide enough systems 
was indeed a strong motivation to considering under-resolved DNS as a modeling strategy 

m 

3.2 From temporal to spatiotemporal dynamics through prob- 
abihty distributions 

From a statistical point of view, energy fluctuations are best characterized by their proba- 
bility distribution functions (PDFs). Figure [8] plots, in lin-log coordinates, the normalized 
histograms Il{E) of E{t) defined as the volume average of |v^, i.e. E = |(t/rms)^ and 
recorded for four different domains, L = 16.9, 35.3, 76.2, and 98.5 at different values of 
i?, during runs of duration up to 3 • 10"^ time units. This is much longer than for Fig. [5] 
or [71 and sufficient for the present purpose since the remark about the duration of the 
simulation does not apply as long as the system is not on the verge of decaying. All the 
curves display a marked hump with some variations. For the smallest domain, FiglHtt, 
both the mean value (see also Fig. [5]) and the most probable value (MPV) of E{t) do 
not change significantly with R] the tail present at large E is presumably a signature 
of the underlying chaotic dynamics. For L = 35.3 in Fig. [8b the curves are roughly 
parabolic (which corresponds to Gaussian distributions), except for R = 380. The MPV 
stays fixed with R like for L = 16.9 whereas the extended tail towards smaller E for 
R = 380 suggests approaches to the laminar state which, when sufficiently marked, lead 
to turbulence breakdown, so that the system likely not far above the value below which 
decay can happen in a short time. The same tendency is observed in Fig. [Ht for an larger 
domain, L = 76.2, but now a downward shift of the MPV of E is observed, which is best 
attributed to the existence and growth of the laminar fraction as R is decreased, laminar 
fraction that will be more conspicuous at larger L as shown in Fig. [71 The PDF also 
displays a low-end tail for R = 360 which again implies that turbulence breakdown can 
take place if one is patient enough in pursuing the simulation, and will do so more easily 
at a lower R. The last case is for L = 98.5 in Fig. [8)1, showing cleaner parabolic shapes 
and the systematic shift in the MPV of E. Since R = 330 is sufficiently above Rg = 325, 
the system is not at risk to decay and the PDF has no low-E" tail. 

From these curves, one main feature emerges: the position of the MPV may be variable 
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Figure 8: (Color online) Probability distribution functions Il{E) as functions of E for 
different system sizes; values of R are indicated in the legends, (a) L = 16.9; (b) L = 35.3; 
(c) L = 76.2; (d) L = 98.5. Notice the different scales for E on the horizontal axis. 
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Figure 9: Width a of the PDFs of E compensated for the presumed extensive behavior 
of the featureless turbulent regime at i? = 420 for L = 16.9, 35.3, 76.2, and 98.5, from 
data presented in Fig. [8l 

or not. Snapshots of the solutions show that the decrease of the MPV with R reflects the 
presence of a larger and larger region of the system which has returned to laminar flow, 
temporarily or persistently, and the PDFs tell us that it is in a statistically signiflcant 
way. We interpret this as the signature of the change in the dynamics from temporal to 
spatiotemporal and conclude that one can decide this issue blindly from the consideration 
of time series of E only. From the variation of the MPVs, one can locate the change for 
L around 76.2. A secondary feature is the presence or not of low-end tail in the PDF of 
which is an indication of the robustness of the state considered. In the absence of a 
tail, the distribution is essentially Gaussian and the system is not at a risk to decay at 
the given value of R. If an exponential tail is present, it means that the decay probability 
is small but signiflcant, which means that if the experiment is long enough, the flow will 
decay. 

Let us now consider the featureless turbulent regime at i? = 420 where the PDFs 
displayed in Fig. [8] are all nicely one-humped. A quantitative characterization of their 
shapes is obtained by fltting lnn(£') against a polynomial in the form aE'^ + bE + c. 
The width of the distribution can be deflned from a as a := |a|~^/^. By eye, it can be 
seen in Fig. [8] that a decreases as L increases, which can be understood as a consequence 
of averaging over wider domains. Assuming local random behavior at the scale of the 
MFU, one can view the system as an assembly of independent sub-systems (which is of 
course not the case but may serve as a template), in which case one expects the standard 
deviation of the fluctuations to vary as the inverse square-root of their number, i.e. as 
the inverse of the square-root of its surface. Accordingly, in Fig. [9] we plot a compensated 
by its presumed variation with the size, i.e. a^LxL^ as a function of v^L^^- The fact 
this quantity is indeed approximately constant supports the underlying assumption but, 
by virtue of contrast, also suggests that it will have to be reexamined when the pattern 
will set in as a result of the interaction between MFUs at lower R. 

3.3 Correlation lengths and the temporal-spatiotemporal issue 

In order to gain further insight into the temporal spatiotemporal transition at a 'mi- 
croscopic' level, we now study the changes in the correlation functions upon domain size 
variation. Here we focus on the streamwise perturbation velocity component obtained 
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Figure 10: (Color online) Time- averaged streamwise velocity correlations at i? = 410 for 
various domain sizes as indicated in the legends. (a,c) Along the direction z. (b,d) Along 
the streamwise direction x. Averaging over 1000 successive snapshots taken every 5t = 1 
is performed. 



upon subtraction of the base flow proflle, u = u — y. We consider the streamwise and 
spanwise time-averaged correlation functions in the mid-plane y = where the base flow 
is strictly zero. They are deflned as: 

^ / (72(x + xx,t)72(x,t)) \ r = f + (A) 

where x is the in-plane running point while ( ) and the overlines denote averaging over 
space and time, respectively. Owing to the periodic boundary conditions in z and 
the correlation functions are periodic with periods Lx and L^. Furthermore, symmetry 
considerations lead to expect them to be symmetric with respect to the origin, which 
is approximately the case. Accordingly, we only display one-sided symmetry-averaged 
correlations proflles over [0, 1/^/2] or [0, Lx/2] for the different domains that we have 
considered. They are presented for i? = 410 in Figure [TOl which shows and Cx^ in 
panels (a,c) and (b,d), respectively. 

Due to the constraints brought by the boundary conditions, both and Cx show that 
small domains are highly correlated. It should however be noticed that Cz passes zero 
even for the smallest domain, which is consistent with the presence of a pair of streaks 
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(see Fig. [3j). Information about processes at the scale of streak, such as the regeneration 
mechanism of turbulence can thus be studied in such domains. That the process depends 
much more on the velocity profile variation in z than in x can be understood from the 
high level of correlation along x (Fig. [TOb . L=16.9). With increasing domain size L, 
Cz oscillates with a period typical of the streaks. Oscillations seem to persist even at 
L = 143.1. On the other hand, the streamwise correlation function decreases much 
more slowly and may reach zero, testifying for the very elongated character of these 
structures. 

The above discussion shows that L = 76.2-87.3 forms some kind of a boundary be- 
tween the two distinct regions which exhibit temporal and spatiotemporal dynamics. This 
boundary will be further explored in section 13.41 below. It is suggested that only domains 
able to accommodate these naturally existing elongated streaks are those that will exhibit 
banded patterns at lower R. In this connection let us notice that, in their work, Barkley 
and Tuckerman were implicitly taking this feature into account by using their tilted do- 
mains since periodic conditions were correlating streaks that were shifted by one spanwise 
MFU width Xz so that their streamwise amplitude could be modulated on a scale much 
longer than the period imposed by their pseudo-streamwise conditions at L^^/ = A^/ sin0, 
where 9 is the angle between their domain short direction and the x axis [2j . In the simu- 
lation corresponding to their Fig. 1, modulation took place over 9 periods, corresponding 
to simulations in a periodic box of size Lx = 97 and = 43 aligned with the flow. The 
limitation of their approach is that correlation in the spanwise direction is somewhat en- 
hanced by the shift condition, which ends in patterns that look much more regular than 
what one obtains in boxes like our largest ones, L = 109 {Lx = 100, = 45, close to 
their effective dimensions) or L = 143.1 (see snapshots in Fig. [71). 

3.4 Bifurcation diagram of transition in plane Couette flow 

Bifurcation diagrams for PCF at different values of L are displayed in Fig. [HI left panel. 
In the right panel we show them unfolded according to the system size, i.e. in a R-L-D 
coordinate system, the presentation adopted for Fig. [3 The diagrams are all obtained 
following the same adiabatic protocol with AR = 10 and AT = 10^, starting from states 
prepared at i? = 420, except for systems with L = 16.9 and 25.6 that are started at 
R = 450. The distance D used to characterize the state of the system is t/rms averaged 
over the last 600 time units. In all cases, an abrupt transition from sustained turbulent 
regime to laminar flow is observed as R is decreased. It should however be stressed that 
the exact quantitative position of the jump depends on the protocol: 

• When L is small, the transition is probabilistic [iTH [35j but a single trajectory 
is followed. A better determination would necessitate to determine PDFs of transient 
lifetimes and the value of R at which the mean lifetime is larger than some value beyond 
which the flow regime would be considered as sustained. 

• At larger L, the relaxation of turbulence has also probabilistic features and choosing 
AT too small may not ensure us that the system has explored a sufficiently large part 
of its accessible phase space so that the result is still sensitive to the initial state at the 
value of R considered, which is the state at the end of the simulation at i? + AR. On the 
other hand, when AR is too large, setting the system at R from R + AR can be a large 
perturbation which may place the initial condition outside the attraction basin of the 
turbulent state at the considered i?, leading to an overestimation of the R corresponding 
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Figure 11: (Color online) (a) Bifurcation diagrams in the R D plane for the different 
system sizes indicated in the legend, with D = t/rms- (b) Three-dimensional unfolding of 
(a) according to system size L. 

to the jump in the bifurcation diagram for that L. This case corresponds to the quench 
experiments in which the initial state is systematically taken to be fully turbulent j6]. 
The risk is more limited if AR is small enough as in the adiabatic protocol that we are 
following. 

We are confident that the decrease of the Reynolds number at which decay to laminar 
flow occurs seen in Fig. [TT]as L increases is not an artifact of the protocol. This decrease 
in R means that, when L is increased, the flow needs to be forced less vigorously to remain 
turbulent because its effective phase space has enlarged so that it has more freedom to 
evolve non-trivially in a spatiotemporal manner, rather than strictly temporally when it is 
conflned by the lateral boundary conditions at small distances. The increase in the region 
of non-laminar flow with increasing domain size can also be inferred from the plot in 
Fig. [T2l where colors from blue = zero to red shows [/rms by increasing values. Also shown 
in this flgure are black-and-white (BW) coded snapshots of the system for a series of 
sizes and Reynolds numbers. The corresponding cases are indicated by white dots in the 
left panel. These images are obtained by coarse-graining ||v|p from the flnal state of each 
simulation at the corresponding R in cells of size 1x1x1, following a procedure described 
in [33]. The top {y > 0) and bottom {y < 0) cells are separately coarse-grained to take the 
lateral shifts of the laminar-turbulent regions in both halves into account |9l [3l [33] . B W 
thresholding is then performed, a pixel being termed 'turbulent' and B-colored when the 
mean of the energy of the two corresponding stacked cells is larger than half the average 
energy in the whole system, 'laminar' and W-colored otherwise. This automated cut-off 
criterion has been found suitable all across the ranges of R and domain sizes studied, 
yielding pictures visually similar to their ^-averaged, non-coarse-grained counterparts. 

The observation of the bifurcation diagrams in Fig. [TT] (left) does not lend itself to 
clearly differentiate the regions where well-oriented patterns occur and where they do not, 
though for domains L = 65.4 and smaller, [/rms is seen not to vary much before decay. 
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Figure 12: Left: (Color online) Bifurcation diagram in the R L plane with color coded 
values of D = t/rms; red (rendered medium gray) correspond to the largest value of D 
(turbulent at large L, chaotic at small L), blue (rendered dark gray) is for laminar flow 
where D = 0^ patterns and rigged states have intermediate values of D from yellow to 
orange (rendered light to medium gray). Above and on the right: BW representation of 
turbulent-laminar patches for L = 143.1 (top), 98.5, 76.2 and 35.3 (bottom). Correspond- 
ing Reynolds numbers are indicated at the bottom right corner of each image. White dots 
in the left panel show the location of the states corresponding to the BW plots. Cases 
(L = 143.1, R = 300) and (L = 98.5, R = 320) are decaying and eventually become 
laminar. 
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Figure 13: (Color online) Streamwise velocity fluctuation u in the plane y = for diflFerent 
domains at i? = 350: L = 76.2, {9 = 23.2°, top), 87.3 {9 = 23.6°, center-left), 98.5 
{9 = 24.0°, center-right), 109.6 {9 = 24.2°, bottom-left), 143.1 {9 = 26.6°, bottom-right). 



whereas in larger domains a gradual change is observed. However, the BW energy plots 
in Fig. [H e.g. (L = 143.1, i? = 400) or (L = 98.5, i? = 380), show that apart from 
the patterned region, there exists a parameter region where the surface of the system is 
'moth-eaten' or 'riddled' with fluctuating laminar troughs. This can also be realized from 
the changing magnitude of [/rms in the color plot if its slow decrease with R is attributed 
to a concomitant increase of the laminar fraction. A hatched domain limited by a white 
dashed line in the left part of Fig. [12] roughly indicates the presence of these riddled 
states. This particular regime thus appears as a precursor to the banded regime when 
the patterns are approached by increasing L or decreasing i?, situated in the bifurcation 
diagram like the intermittent regime described by Prigent et al [32] and Barkley & 
Tuckerman [2j and to which it corresponds closely. A possible equivalent of it in pipe flow 
seems to be the intermittent laminar-flash regime described in f28] near the transition to 
uniform turbulence that we call featureless. In our numerical conflguration, the bifurcation 
diagram reaches its eventual large aspect-ratio limit 'featureless-riddled-banded-laminar' 
for systems of sizes beyond L = 76.2. 

Further visualization in Fig. [T3l shows systems with L > 76.2 where one can visually 
observe the tendency to form a well-organised pattern as L is increased. Below this 
domain size, i.e. L < 65.4, no pattern or troughs could be observed upon reducing R 
down to the final laminar state. Accordingly, patterns seems to need systems larger than 
some minimal size to exist. However, this size should rather be considered to mark some 
crossover than to be 'critical' in the usual sense because, there cannot be any decisive 
criterion to decide whether one has more or less oblique, elongated, laminar troughs like 
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for L = 76.2, or already well-formed bands, like for L = 87.3. As noted above, these 
laminar troughs are also found in very large domains at higher R and are precursors to 
the bands when, at given L, R is reduced adiabatically from well above Rt. In fact, 
from the plots in Fig. [TT] (left), one can see that R = 350 is close to the transition to 
laminar flow for L = 76.2 and 87.3, hence at the corner of the white hatched region in 
Fig. [121 Systems in the crossover region thus go from turbulent to laminar by skipping 
the patterned stage but still going through the intermittent-riddled regime. 

4 Summary and Conclusions 

The formation of banded patterns in flows such as plane Couette flow, circular Couette 
flow and plane Poiseuille flow is typically a spatiotemporal problem. And true enough, 
this phenomenon has been experimentally observed in domains with lateral dimensions 
that are more than two orders of magnitude larger than the wall-normal distance. Numer- 
ical simulations of the flow fleld in these large systems have become possible only recently, 
thanks to the less stringent resolution requirements to resolve all the scales of low-i? tur- 
bulent flow and to our ever increasing computational capabilities. Despite these favorable 
circumstances, simulations cannot be performed over long enough durations to bring def- 
inite answers about the statistics of the transitional regime in the large aspect ratio limit. 
Previously, with fully developed wall-bounded turbulence in mind (high-i?), numerical 
computations and theory were more focused on small domains. In the moderate- i? range 
of interest to the problem of transition, these small systems are best analyzed within 
the framework of flnite-dimensional dynamical systems theory. The present work stays 
in between these two types of studies, wherein we start with low-i? turbulence in small 
domains where temporal dynamics takes place, and increase the domain size gradually to 
reach the spatiotemporal regime. An adiabatic decrease of R for each domain then gives 
information about the different regimes visited by the system from featureless turbulence 
to laminar flow. 

Upon decreasing i?, starting from a featureless fully turbulent flow, we have obtained 
that, in small systems, the mean value of the perturbation energy remains roughly con- 
stant as a function of R down to the point where the chaotic regime breaks down into 
laminar flow, whereas it is regularly decreasing in larger systems. Furthermore, visu- 
alization of the flow fleld indicates that the mean-energy variation corresponds to the 
system entering a spatiotemporal regime that at flrst (high R) presents fluctuating lam- 
inar troughs and later (low R) steady oblique patterns of bands alternately laminar and 
turbulent, but absent in the small systems in which spatial coherence implies a mostly 
temporal dynamics. 

In the featureless regime, fluctuations of the mean perturbation energy have distribu- 
tion probabilities with a single marked hump around the mean for all domain sizes. Their 
variance is decreasing with increasing system size, and fluctuations in the widest domains 
are close to normal, which can be understood as a result of the additivity of local (MFU 
scale) fluctuations that remain little-correlated as long as the system is sufliciently far 
from the laminar breakdown, expressing the extensive character of the featureless regime. 

A study of the spanwise and streamwise dependence of the streamwise velocity com- 
ponent correlation function as functions of the system size has also been performed. The 
spanwise variation accounts for the self-sustaining process by displaying a periodic de- 
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pendence at the scale expected for the streaks. More importantly, the slow streamwise 
variation expresses correlations reinforced by the periodic boundary conditions in the 
smallest systems. The ability to take into account this streamwise dependance properly 
seems an important ingredient in order to explain the absence of patterning noticed in 
these systems. Indirectly, this also explain the observation of the pattern in the Barkley- 
Tuckerman tilted but short domains since their domain width was chosen to fulfill a 
commensurability condition ensuring the periodic continuation of streaks at the tilt angle 
observed in the experiments, thus mimicking longer domains in the streamwise direction. 
In our computations, patterns were found for L > 76.2, i.e. Lx > 70. 

Inspection of fiow fields (e.g. Fig. [13] or the BW insets in Figure [T2l) indicates that, in 
accordance with the information gained from correlation functions, one can only observe 
intermittent laminar troughs for L = 76.2, and see any organized pattern of laminar- 
turbulent oblique bands only for L = 87.3 and beyond. 

In Figure El a conjecture was presented about the bifurcation diagram of PCF, con- 
necting systems of MFU size exhibiting temporal dynamics [3^ to large aspect-ratios 
systems displaying spatiotemporal dynamics and patterns [32]. Substantiating this con- 
jecture at a quantitative level. Figures [H] and [12] display the bifurcation diagrams ob- 
tained for various system sizes. The smallest domains follow the direct route '[chaotic 
flow] [laminar flow]' {via chaotic transients), whereas in very large domains one has 
'[turbulent flow] [riddled regime] [oblique pattern] [laminar flow]'. For L in the 
range 75-85 the [oblique pattern] stage is skipped, which marks the crossover from tempo- 
ral to spatiotemporal dynamics. Furthermore, the threshold at which the pattern decays, 
i?g, is shown to decreases with increasing system size, seemingly tending to a constant in 
the large aspect-ratio limit. A precise quantitative estimate of this limit was however out- 
side the scope of this paper since the flnal turbulence breakdown still keeps probabilistic 
features that requires statistics (longer time series, large number of independent realiza- 
tions), especially in what regards the occurrence of turbulent patches issued fragmented 
bands turning into turbulent spots close to i?g. At any rate the exact value is of little in- 
terest, all the more that periodic boundary conditions tend to stabilize the pattern, which 
thus artiflcially decreases i?g at moderate aspect ratios. The observed trends (decrease of 
Rg as L is increases, shrinking of the band regime to the beneflt of the intermittent riddled 
regime) however go in the same direction as the experimental flndings showing that at 
intermediate aspect ratio, the transitional range is pushed at higher R and that patterns 
can hardly be observed — see [38j and Bottin's thesis [5^ §4.4] — though no quantitative 
link can be made in view of the differences in the lateral boundary conditions. 

All in all, the recognition of a crossover size beyond which the transition to/from 
turbulence in PCF is undoubtedly a spatiotemporal process, and the need of a faithful 
account of streamwise correlations at this scale, either directly or indirectly via the tilted- 
domain trick, seem our most important observations. They might help us to unravel the 
physical mechanism behind the organization of low-R turbulence in wall-bounded flows, 
which still largely remains an enigma. 
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